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ABSTRACT 

Standard stellar evolution theory is inconsistent with the observed isotopic 
carbon ratio, 12 C/ 13 C, in evolved stars. This theory is also inconsistent with 
the 3 He/H abundance ratios observed in Galactic H n regions, when combined 
with chemical evolution theory. These discrepancies have been attributed to an 
extra, non-standard mixing which further processes material during the RGB and 
should lower both the 12 C/ 13 C and 3 He/H abundance ratios for stars with masses 
< 2 M Q . Measurements of isotopic ratios in planetary nebulae probe material 
which escapes the star to be further processed by future generations of stars. 

We have measured the carbon isotopic abundance ratio, 12 C/ 13 C, in 11 plan- 
etary nebulae (PNe) by observing the J=2— >1 and J=3^2 millimeter transitions 
of 12 CO and 13 CO in molecular clouds associated with the PNe. A large velocity 
gradient (LVG) model has been used to determine the physical conditions for 
each PNe where both transitions have been detected. We detect both 12 CO and 
13 C0 in 9 PNe. If 12 CO/ 13 CO = 12 C/ 13 C, the range of 12 C/ 13 C is 2.2-31. Our 
results support theories which include some form of extra mixing. 

Subject headings: ISM: abundances - - ISM: molecules - - planetary nebulae: 
general — radio lines: ISM 



Introduction 



Measuring isotopic abundance ratios is one method of probing how stars of various 
masses process material and how the Galaxy redistributes this material over time — that is, 
stellar and Galactic chemical evolution. Some of the light elements H, 2 H, 3 He, 4 He, and 
7 Li, were produced during the era of primordial nucleosynthesis, so these ratios provide 
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additional constraints for Big Bang nucleosynthesis. H n regions or molecular clouds are 
examples of zero-age objects which are young relative to the age of the Galaxy. Therefore 
these abundance ratios chronicle the results of billions of years of Galactic chemical evolution. 
Observations of isotopic ratios in stars and planetary nebulae (PNe) measure locally how 
stars influence chemical evolution. In particular, PNe probe material which has been ejected 
from low-mass (M < 2M ) to intermediate-mass (M ~ 2 — 5M ) stars to be further 
processed by future generations. 

Measurements of 3 He in the Galaxy have led to what has been called "The 3 He Problem" 
(Galli et al. 1997). Standard stellar evolution theory predicts that 3 He is produced in 
significant quantities by stars of 1 — 2 M and that PNe with progenitors in that mass range 
should have 3 He/H ~ several xl(T 4 (Iben 1967; Rood 1972; Vassiliadis & Wood 1993; 
Weiss et al. 1996). (All abundance ratios will be by number unless otherwise noted.) 
Galactic chemical evolution models which integrate over a range of initial masses predict 
that 3 He is produced in prodigious amounts (Rood et al. 1976; Dearborn et al. 1986). A 
survey of 6 PNe produced a detection of 3 He in NGC3242 (Rood, Bania, & Wilson 1992; 
Balser, Rood, & Bania 1999) and probable detections in two other sources (Balser et al. 
1997) with 3 He/H ~ 10~ 4 - 10 -3 as predicted by standard stellar evolution models. This 
PNe sample was highly biased, however, and was selected to maximize the likelihood of 
a detection. Observations of 3 He in protosolar material (Geiss 1993), the local interstellar 
medium (Gloeckler & Geiss 1996), and Galactic H n regions (Balser et al. 1999) all indicate 
a 3 He/H abundance ratio of ~ 2 x 10 -5 . Although H n regions are expected to be zero-age 
objects there is no evidence for any stellar 3 He enrichment during the last 4.5 Gyr (Bania et 
al. 2002). The observed abundances of 3 He throughout the Galaxy when combined with 
the observations that some stars do indeed produce 3 He strongly disagree with standard 
chemical evolution models (Steigman & Tosi 1992; Vangioni-Flam et al. 1994; Steigman & 
Tosi 1995; Dearborn et al. 1996; Prantzos 1996; Fields 1996; Scully et al. 1996; Galli et 
al. 1997). 

Anomalies in the 12 C/ 13 C abundance ratio in red giant stars has stimulated studies of 
non-standard stellar evolution models (Charbonnel 1994, 1995; Wasserburg et al. 1995; 
Weiss et al. 1996; Sackmann & Boothroyd 1999; Boothroyd & Sackmann 1999). All of 
these models include some form of extra mixing which transports material in the convective 
envelope to lower regions which are hot enough for further nuclear processing and then 
transports it back up into the convective envelope. The extra processing that occurs is 
called cool bottom processing (CBP). 

During the first dredge- up on the red giant branch (RGB), the convective envelope 
reaches regions where 12 C had previously been processed into 13 C and 14 N; thus we expect 
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12 C/ 13 C to decrease from ~ 90 to 25 for both low-mass and intermediate-mass stars (Iben 
1964; Iben & Renzini 1984). Observations of carbon in the surfaces of red giant stars reveal 
lower 12 C/ 13 C ratios (5-20) for stars below 2M than theoretically expected (e.g., Gilroy 
1989). Studies which measure the evolution of the surface 12 C/ 13 C as a star ascends the RGB 
in M67 (Gilroy & Brown 1991) and in field stars (Charbonnel et al. 1998) indicate that 
additional processing occurs between the base and tip of the RGB. This has been attributed 
to extra mixing below the convective envelope. 

During the asymptotic giant branch (AGB) the theoretical abundances depend on the 
interplay between dredge-up events, mass loss, hot bottom burning, and non-standard mixing 
(CBP). During the third dredge-up on the AGB significant pollution of 12 C will increase the 
12 C/ 13 C ratio if hot bottom burning does not take place (Forestini & Charbonnel 1997). 
In low-mass stars, 12 C/ 13 C may decrease to ~ 4 if AGB cool bottom processing occurs 
(Wasserburg et al. 1995). For AGB stars more massive than ~ 4M , hot bottom burning 
drives the 12 C/ 13 C towards its equilibrium value of 3.5 (Frost et al. 1998). 

Observations of abundances in PNe measure the cumulative effects of different mixing 
and nuclear processing events during the entire evolution of their progenitors. In general 
the abundance of 3 He should follow the 12 C/ 13 C ratio; that is, if 3 He is destroyed then 
the 12 C/ 13 C ratio should decrease. Therefore given the current observational challenges of 
detecting 3 He in PNe, an alternative approach is to measure 12 C/ 13 C. 

There have been several studies of 12 C/ 13 C in the molecular gas in PNe using millimeter 
transitions of CO (Bachiller et al. 1989; Jaminet et al. 1991; Cox et al. 1992; Bachiller et 
al. 1997; Palla et al. 2000). In these studies it is assumed that the 12 C/ 13 C ratios measured 
in molecular clouds near the PN refer to the abundances in the PN itself. Alternatively, 
Clegg et al. (1997) and Miskey et al. (2000) have used the 13 C 3 1/2 P° - } /2 S transition at 
1909.597A to measure the 12 C/ 13 C ratio in several PNe. From these studies the measured 
12 C/ 13 C abundance ratios are between ~ 10 to 40. Palla et al. (2000; hereafter PBSTG) 
concluded that in the majority of PNe the 12 C/ 13 C ratio is lower than expected from standard 
stellar models and are consistent with non-standard mixing models which incorporate CBP. 

In this paper we present new observations of 12 CO and 13 CO which enlarge the PN 
sample. Usually isotope ratios are estimated using optically thin lines, however, since the 
CO intensities are small, and since the 16 0/ 18 ratios may vary, the 12 CO and 13 CO lines 
are used to deduce the 12 C/ 13 C ratio. Because of the resulting uncertainty in opacity, it 
is important that at least two transitions of CO be observed to constrain radiative models 
used to calculate the 12 C/ 13 C abundance ratio. The observations are discussed in §2 and 
the general results are in §3. In §4 the 12 C/ 13 C abundance ratios are calculated using the 
assumption that the optical depths of both 13 CO and 12 CO are small and that the transitions 
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are themalized. More sophisticated methods are explored in §5. In §6 the implications with 
respect to stellar and Galactic evolution are discussed. A summary is given in §7. 

2. Observations 

The observed sample of PNe are listed in Table 1. The PNe were selected from the CO 
survey of Huggins & Healy (1989) based on the estimated intensity of 13 CO for an expected 
12 CO/ 13 CO ratio of 25. Listed are the source name, the Galactic coordinates rounded to the 
nearest degree, the equatorial coordinates, the observed offset position, the LSR velocity, 
the angular size, the progenitor main sequence mass, and the morphology. The progenitor 
masses, M ms , are taken from PBSTG for objects in common or determined by first using the 
central star masses determined by Gorny, Stasihska, & Tylenda (1997) or Stasihska, Gorny, 
& Tylenda (1997) and the initial-final mass relation of Herwig (1996). The PNe morphology 
is defined as round (R), elliptical (E), or bipolar (B) (Balick 1987; Manchado et al. 1996; 
Huggins et al. 1996). 

An observing log is summarized in Table 2 which includes the spectral transitions, the 
telescope's half-power beam-width, and the velocity resolution. Observations of the 12 CO 
and 13 CO (J=2— >1) rotational transitions were made with the National Radio Astronomy 
Observatory 5 12 m telescope on 1997 May and 1997 December. The 200-300 GHz dual- 
polarization SIS receiver was employed. Both the filter-bank and hybrid spectrometers were 
used simultaneously and were configured to provide spectral resolutions of 1 MHz and 2 MHz 
per channel which produce a velocity resolution of 1.4 km sec -1 and 2.7 km sec -1 per channel, 
respectively. The spectra were obtained by beam switching plus position switching (BSP). 
The reference position (OFF) was typically 4 arcmin from the source position (ON). The 
pointing should be accurate to within ~ 5 arcsec. The data were calibrated on-line to the 
radiation temperature (T^) scale which corrects for atmospheric attenuation, radiative loss, 
and rearward and forward scattering and spillover (Kutner & Ulich 1981). In this paper 
we report all intensities in units of the main beam brightness temperature, T mb , using an 
efficiency of rj^ = 0.5 which converts the radiation temperature to the main beam brightness 
temperature (Mangum 1999). 

Observations of the 12 CO and 13 CO (J=3— >2) transitions were made with the Hen- 
rich Hertz Telescope (HHT) 6 during February 1998 and May 1998. The facility 345 GHz 
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dual-polarization SIS receiver was employed, although because of technical difficulties only 
one channel was used. The spectrometers consisted of a 1 GHz, 1024-channel AOS and 
a 250 MHz, 1024-channel AOS with spectral resolutions of 1.0 MHz (0.87kmsec _1 ) and 
0.40 MHz (0.35 km sec -1 ) per channel, respectively. Beam switching was used to produce 
the spectra using a procedure called WSWITCH with a reference offset from the source by 4 
arcmin. (This procedure is effectively identical to BSP at the 12 m.) The pointing model was 
determined to be accurate to within 1-2 arcsec for most locations in the sky. The data were 
calibrated on-line to the T^ temperature scale which corrects for atmospheric attenuation. 
The line intensities were scaled to main beam brightness temperatures using several standard 
calibration sources (x Cyg, R Aql, and IRC+10216). These intensities are consistent with 
an assumed efficiency of r] m b = 0.55, if the mixer sideband ratios are equal. We estimate 
that the uncertainties in the intensity scale are between 15-20% (Wilson et al. 2001). 



3. Results 

The results are summarized in Table 3. Listed are the source name, the transition, the 
peak intensity (T pea k) , the root-mean-squared (RMS) variations in the spectral baseline, the 
velocity of the centroid (V ccntr oid)) the full-width half-maximum of the linewidth based on a 
single-component Gaussian fit (AV), and the integrated intensity (/). Only AV^ is based on 
a Gaussian model. The data were analyzed independently using the spectral line analysis 
software CLASS, UniPops, and COMB. The majority of the sources were detected in both 
the 12 CO and 13 CO J=2 — >1 transitions. Four objects were also detected in the J=3^2 
transitions which require higher excitation. 

Figure 1 summarizes the spectra taken with both telescopes. The first two columns 
correspond to the NRAO 12 m J=2^1 transitions while the last two columns correspond to 
the HHT 10 m J=3^2 transitions. The spectra have been shifted in velocity to align the 
profiles relative to the LSR velocity in Table 1. For several sources there is evidence for 
complicated velocity structure (e.g. multiple peaks, asymmetries, etc.). In these cases this 
structure is being averaged over the telescope's beam and thus our analysis refers to the bulk 
properties of the molecular emitting gas. 

For PNe observed at both transitions we measure the J=3^2 transition to be signifi- 
cantly more intense than the J=2— >1 transition. This effect is expected based on solely the 
radiative transfer in gas with physical conditions characteristic of PNe. 
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In Figure 2, the radiation temperature is plotted versus the upper J level rotational 
transition for CO over a range of densities. Across a broad range of physical conditions, the 
J=3^2 transition produces the brightest intensity; only in the lowest volume density-column 
density models is the J=2^1 intensity higher. The generally better coupling between smaller 
beams (obtained at higher frequencies) and PN structure also suggests that CO searches in 
PNe may be better carried out in the sub-millimeter regime; atmospheric and hardware 
limitations may mitigate this effect. 



4. 12 C/ 13 C Abundance Ratios 

The derived 12 C/ 13 C ratios for our PN sample are listed in Table 4. In columns 3 and 
4 the 12 C/ 13 C ratios are shown for the J=2^1 and J=3^2 transitions, respectively. Values 
based on the integrated intensities are expected to be more accurate, but values based on 
the peak intensities are also given (in parentheses). The 12 CO/ 13 CO abundance ratios are 
believed to be a direct measure of the 12 C/ 13 C ratio; therefore we have assumed that 12 C/ 13 C 
= 12 CO/ 13 CO. Chemical fractionation will lower the observed 12 CO/ 13 CO ratios, but this 
effect should be reduced by the observed high temperatures in the molecular gas (~ 25 — 60 
K, Bachiller et al. 1997). Selective photo-dissociation of 13 CO causes an increase in the 
12 CO/ 13 CO ratio but this will be counterbalanced by the charge exchange reaction 13 C + + 
i2 C0 i2c+ + ^CO (Likkel et al. 1988). 

The range of 12 C/ 13 C ratios are lower than the values predicted by standard stellar 
models. About half of the objects observed in this study were also included in the sample of 
PBSTG. In most cases our ratios agree with the PBSTG values shown in the last column in 
Table 4. The peak main brightness temperatures for our observations are lower in all cases, 
as expected for the larger beam in our study. However the intensity is not proportional to 
the beam areas. 

These 12 C/ 13 C ratios were calculated under the assumptions that the filling fraction for 
the gas is the same for both isotopomers, both rotational transitions are thermalized to the 
same temperature, and that the emission lines are optically thin. The first assumption seems 
reasonable based on the spectra shown in Figure 1. A comparison of the different transitions 
reveals nearly identical Vlsr and line shape characteristics, indicating that similar volumes of 
gas are sampled by the two transitions. The excitation characteristics of the low-J transitions 
of CO also supports the assumption of a thermalized, uniform temperature. The small dipole 
moment of CO results in a low critical density for rotational excitation; thus most PNe have 
densities high enough for significant CO excitation. The effect of line opacity in CO may be 
accounted for in multi-transition models of radiation transport. 
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5. Large Velocity Gradient Models 

In order to better interpret our results, the radiative transfer in the PNe is examined 
through the use of the Sobelev or Large Velocity Gradient (LVG) approximation (Sobolev 
1958, 1960). This model provides a numerical solution to the coupled equations of radiative 
transfer and statistical equilibrium for the rotational transitions of CO (see e.g. Leung & 
Liszt 1976). In particular, we used the RAD package which is a part of the MIRIAD 7 software 
package (Mundy 1992). Only the four PNe for which we have measured multiple transitions 
of 12 CO and 13 CO emission are considered. For each source a two-dimensional grid composed 
of volume density (n) and column density (N), over a range of physical conditions relevant to 
molecular gas in PNe, is generated. We adopt a kinetic temperature = 60 K. This value 
is towards the high end of temperatures suggested by Bachiller et al. (1997); however, the 
higher J transitions of CO typically trace warmer regions and the LVG grids are consistent 
with T k = 60 K. Data taken with the HHT (J=3^2) have been spatially smoothed to match 
the resolution of the 12 m (J=2^1) data. We assume that the emission fills the beam for 
both pairs of transitions, or has a similar beam filling fraction. We obtained y 2 values from 
a comparison of the data to the model grids as a measure of the quality of the fits. 

The derived 12 C/ 13 C ratios are listed in column 5 of Table 4. The quoted uncertainties 
were determined by propagating the measured error in the integrated intensities from Table 
3. These statistical uncertainties may be much smaller than the systematic uncertainties. 
For the PNe for which there only exists J=2— >1 detections the integrated intensity ratios are 
adopted. Where the line shape is parabolic the line is optically thick so the 12 C/ 13 C ratios 
are cited as lower limits. For PNe where we have detected both the J=2^1 and J=3^2 
transitions an LVG analysis is performed and a 12 C/ 13 C ratio determined. Below we discuss 
each of these sources in more detail. 

AFGL 618 — Using 13 CO we derive n = 6.8 x 10 4 cm _3 and N = 6.8 x 10 15 cm~ 2 . Using 
these values, we obtain N = 3.2 x 10 16 cm -2 for 12 CO, yielding a ratio of 12 CO/ 13 CO ~ 5. 
Overall this value is consistent with the observations. This model overproduces the J=2^1 
emission by about 20% and the best fit density to the J=2— »1 is a factor of ~2 lower than that 
derived from the 13 CO. This may suggest that the observed 12 CO/ 13 CO ratio is a lower limit, 
since the 12 CO and 13 CO emission regions are not identical, or that density, temperature, 
and abundance gradients exist which are not isolated in velocity space. Although the derived 
opacity in the J=3^2 and J=2^1 lines from the LVG grids is low, the line shape is strongly 
suggestive of high opacity. We therefore cautiously cite the adopted ratio as a lower limit 
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only. 

Ml-16— The model for Ml-16 produces n = 6.8 x 10 4 cm- 3 and JV = 6.8 x 10 15 cm- 2 
for 13 CO, and JV = 1.5 x 10 16 cm~ 2 for 12 CO, yielding 12 CO/ 13 CO = 2.2, very similar to that 
derived from both the peak and integrated intensities. Although this value is very low, the 
model suggests that the opacity remains low for both transitions (e.g., r < 0.05 for 12 CO). 
PBSTG consider this source to be a lower limit, however, presumably because of high opacity. 
Examining the shape of the spectrum, the profile seems most like the classical flat topped 
structure of an optically thin, unresolved shell. Therefore we adopt our model-derived value 
for the ratio in this source. 

IRAS 21282+5050— Using 13 CO we obtain n = 3.2 x 10 5 cm" 3 and TV = 1.0 x 10 15 cm" 2 . 
This produces a column density of 3.2 x 10 16 cm -2 and a ratio of 12 CO/ 13 CO = 32. The 
volume density for IRAS 21282+5050 is larger than that obtained from the other sources, 
perhaps indicative of its youth. The shape of its spectral profile is difficult to discern given 
the signal-noise in our data. For 12 CO in the J=3^2 line the opacity should be higher but 
this line seems more parabolic than flat topped. Moreover, the values obtained from the 
ratios of the integrated intensities and peak intensities are not consistent, most likely due to 
the low signal-to-noise ratio in the 13 CO data. Due to these concerns, we adopt 12 CO/ 13 CO 
= 32 as a lower limit. 

NGC 7293 — Numerous models were run for NGC 7293 without obtaining any reasonable 
fits to the individual 12 CO or 13 CO data. In addition, the derived physical characteristics 
for 13 CO were highly inconsistent with those from 12 CO. Since NGC 7293 is an evolved 
PN, the simplest explanation for these results is that structure in the CO line emitting gas 
couples differently in the two beams, so the J=2^1 and J=3^2 regions are distinct and 
different. The emission profiles from the lines do not suggest a high opacity, however. Since 
the J=3-^2 13 CO transition only has a signal-to-noise ratio of ~ 4 we adopt the J=2 — >1 
ratio. 



6. Discussion 

The carbon isotopic ratio, 12 C/ 13 C, along with other elements are processed during 
different stages of stellar evolution such as the red giant branch (RGB), the asymptotic 
giant branch (AGB), and the planetary nebula (PN). Observations have been made during 
these different stages of evolution which can be used to constrain stellar evolution theory. 

Observations of anomalous 12 C/ 13 C abundance ratios in RGB field stars (Tomkin, Lam- 
bert, & Luck 1975; Tomkin, Luck, & Lambert 1976) and open clusters (Gilroy 1989) 
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have stimulated theoretical studies which introduce an additional mixing process to reduce 
12 C/ 13 C in stars (Sweigart & Mengel 1979; Charbonnel 1994, 1995; Wasserburg et al. 
1995). A detailed study of M67 revealed that the anomalously low 12 C/ 13 C ratios occurred 
only for stars at the tip of the RGB (Gilroy & Brown 1991). Reanalysis of 191 field and 
cluster giants stars using HIPPARCOS parallaxes has determined that 96% of evolved stars, 
objects which have passed the luminosity bump, have 12 C/ 13 C ratios lower than standard 
models predict (Charbonnel & Do Nascimento 1998). Gratton et al. (2000) have measured 
12 C/ 13 C in 62 field metal-poor stars, including 43 additional stars from the literature. They 
determine two distinct mixing episodes. The first dredge-up which produces 12 C/ 13 C ratios 
consistent with standard models ( 12 C/ 13 C = 20-30); and a second mixing episode occurring 
after the luminosity bump, when the molecular weight barrier is removed, where 12 C/ 13 C 
is reduced to 6-10. These results are consistent with the theoretical models which predict 
extra mixing (Charbonnel 1994, 1995). 

Additional processing of carbon and other elements occurs on the AGB. Observations 
have been made of 12 C/ 13 C in AGB stars using mm- wave, optical, and IR spectroscopy. CO 
mm-wave observations of the circumstellar envelops of carbon and oxygen-rich AGB stars 
yield 12 C/ 13 C values between ~ 10-50 (Wannier & Sahai 1987; Kahane & Jura 1994; 
Greaves & Holland 1997). Evidence for some deep mixing process is suggested from an 
optical study of carbon stars where Li-rich objects typically have low (< 15) 12 C/ 13 C ratios 
(Abia & Isern 1997). IR spectroscopy of cool carbon stars produces 12 C/ 13 C ratios between 
30 - 70 with a few very low ratios, 12 C/ 13 C - 4 (Lambert et al. 1986). The high 12 C/ 13 C 
values are expected from third dredge- up of 12 C, while depending on the mass of the star 
the very low 12 C/ 13 C ratios could be explained by hot bottom burning or by deep mixing. 
In contrast Ohnaka & Tsuji (1996) determine 12 C/ 13 C ratios that are 2-3 times smaller 
than Lambert et al. even though there is a significant overlap in the stars observed. There 
have been heated discussions about the analysis of spectra using the iso-intensity method 
(de Laverny & Gustafsson 1998; Ohnaka & Tsuji 1998; de Laverny & Gustafsson 1999; 
Ohnaka, Tsuji & Aoki 2000). Nevertheless, it appears that at least some objects do have 
reduced 12 C/ 13 C ratios compared to that expected from standard models and that some 
extra mixing mechanism is at work. 

Observations of 12 C/ 13 C in PNe measure how carbon has been processed through both 
the RGB and AGB phases. Because this processed material escapes into the interstellar 
medium these 12 C/ 13 C ratios can be used to constrain Galactic chemical evolution models. 
Figure 3 summarizes the current status of 12 C/ 13 C abundance ratios in PNe. Plotted are 
the determined 12 C/ 13 C ratios versus the estimated main sequence progenitor mass (M ms ). 
Only PNe for which progenitor masses could be found are shown. The results from other 
CO observations are taken from PBSTG and displayed as triangles. Measurements using the 
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CIII] transition are shown as stars. Notice that there are five PNe that we have in common 
with PBSTG. These objects are indicated by open triangles. 

For clarity only the (statistical) uncertainties in the 12 C/ 13 C ratio are shown in Figure 
3. The uncertainties in the derived progenitor masses are significant and are mostly due to 
estimates of the PN distance and calculations of the progenitor mass from the current central 
star mass using some initial mass-final mass relationship. These uncertainties can be as large 
as 50% (e.g, Galli et al. 1997; PBSTG). Therefore care must be used when interpreting these 
results. The error bars in Figure 3 for 12 C/ 13 C ratios are based on the uncertainties in the 
measured integrated intensity for each transition. Note that the systematic errors may be 
much larger than these, in some cases. 

Our results are consistent with PBSTG with the exception of NGC6720 where our 
12 C/ 13 C ratio is a factor of two smaller. However, the signal-to-noise ratio for the 13 CO 
J=2^1 transition is only 5 and we do not have multiple transitions for this PN. Therefore 
the 12 CO transition may be optically thick and our 12 C/ 13 C value a lower limit. However 
from the line shapes these data are consistent with optically thin emission. 

The 12 C/ 13 C ratios range between 2.2-31 with progenitor masses from ~ 1 — 4 M . This 
indicates that at least some stars do undergo some non-standard processing such as cool 
bottom processing (CBP). Because of the large uncertainties in determining the progenitor 
mass it is difficult to make any strong statements about the relationship of 12 C/ 13 C with M ms . 
Since CBP is expected to operate only for stars with M ms < 2 M Q the 12 C/ 13 C ratios should 
be larger for stars above this mass. Above 4M hot bottom burning can drive the 12 C/ 13 C 
ratio towards it equilibrium value of 3.5 (Frost et al. 1998). The low 12 C/ 13 C values between 
2-4 M Q might indicate further processing on the AGB or that we have underestimated the 
12 C/ 13 C ratios. It is interesting to note that the lowest 12 C/ 13 C values are classified as 
bipolar PNe. Bipolar PNe tend to lie closer to the Galactic plane with a wider distribution 
of central star masses and have higher N/O ratios (Gorny, Stasihska, & Tylenda 1997) 

Since most of the PNe in our sample have only been detected in one transition we do not 
have sufficient constraints to determine the CO line opacity. If the 12 CO line is optically thick 
then the 12 C/ 13 C values will be underestimated and the points in Figure 3 must be increased. 
Clearly more 12 C/ 13 C ratios are required to improve the statistics, together with multiple CO 
transitions to constrain radiative transfer models. Because of the higher kinetic temperatures 
in PNe, observations in the sub-mm regime may be more useful for this purpose. 

A significant effort has been made to detect CO toward the six PNe in which 3 He have 
been measured to directly establish the relationship between non-standard mixing in low- 
mass stars and the 3 He and 12 C/ 13 C abundances (PBSTG). The stellar models predict that 
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non-standard mixing (CBP) will lower both 3 He/H and 12 C/ 13 C. Based on chemical evolution 
models more than 90% of low-mass stars must undergo CBP to reconcile the observations 
of 3 He in H n regions (Charbonnel et al. 1998). Only one object with a measured 3 He/H 
ratio, NGC6720, has been detected in CO. Since the measured 3 He upper limit of 3 He/H 
= 0.54 ± 0.29 x 10~ 3 is consistent with both the standard model and with models which 
include CBP, the current data do not provide any interesting constraints. Note that the 
progenitor mass of NGC 6720 is slightly above 2 M Q which is where CBP should turn off. 

CO is typically detected in young PNe where the nebula has not evolved sufficiently 
to dissociate the molecular material. In contrast, 3 He is best detected in evolved, older 
PNe which tend to fill the telescope's beam. Therefore it is not a coincidence that there 
is little overlap between the two samples. As pointed out by PBSTG, observations using 
the CIII] transition toward PNe will have a better opportunity to match the 3 He samples. 
Currently, the number of PNe with either CIII] or 3 He detections is sparse. The Green 
Bank Telescope will offer a significant improvement in our ability to detect 3 He + in PNe and 
thus increase the current sample of one well determined ratio (NGC 3242) and two probable 
detections. Observations of CO in PNe will complement any CIII] data. Because of the large 
uncertainties in estimating the progenitor masses we feel that the number of detections must 
be significantly increased to build a strong observational case for non-standard mixing in 
PNe. It is interesting to note that even if CBP occurs in most low-mass stars it will not be 
a significant factor in the chemical evolution of 12 C/ 13 C in the Galaxy since the production 
of carbon is dominated by intermediate-mass stars (PBSTG). Therefore, measurements of 
3 He is one of the few methods to probe how low-mass stars influence the chemical evolution 
of the Galaxy. 



7. Summary 

The following are the main results of the J=2 — >1 and J=3^2 millimeter transitions of 
12 CO and 13 CO in 11 Galactic planetary nebulae (PNe). 

1. We have detected both 12 CO and 13 CO emission in 9 PNe for the J=2^1 transition 
using the NRAO 12 m telescope and in 4 PNe for the J=3— >2 transition using the HHT 10 m 
telescope. Assuming 12 C/ 13 C = 12 CO/ 13 CO the ratios range from 2.2-31. These results 
support an extra mixing process in at least some stars. 

2. Large velocity gradient (LVG) models have been applied to the 4 PNe where both 
transitions have been detected. These models were primarily used to determine if the 12 CO 
line is optically thick, an effect which would overestimate the 12 C/ 13 C ratio. A very low 
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12 C/ 13 C ratio of 2.2 was determined for Ml-16 which does not appear to be a lower limit. 

3. The J=3— »2 transition is significantly more intense than the J=2— »1 transition. This 
result has been confirmed by models using typical densities measured in PNe. This suggests 
that CO observations in PN may be better carried out in the sub-millimeter regime. 

We thank the technical support staff at both the NRAO 12 m telescope and the SMTO. 
In particular, we thank Harold Butner for help during our observations at the HHT. We 
acknowledge fruitful discussions with Corinne Charbonnel and Bob Rood about stellar evo- 
lution. Lastly, we thank Daniele Galli for communicating his results before publication. 
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Fig. 1. — NRAO 12 m spectra of the 12 CO and 13 CO emission in the J=2 — >1 transition 
(first two columns) and SMTO HHT spectra of the 12 CO and 13 CO emission in the J=3^2 
transition (last two columns) for our sample of PNe. The intensity scale has been converted 
into main beam brightness temperature in units of Kelvin. The vertical grey line marks the 
Vlsr of the source. 
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Fig. 2. — Plot of radiation temperature versus lower J transition over a range of physical 
conditions. These LVG models indicate that the J=3^2 line is the peak over the range of 
values relevant to most planetary nebulae. 
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Fig. 3. — The 12 C/ 13 C abundance ratios versus the progenitor main sequence mass. The 
circles are taken from this paper, the triangles are from PBSTG, and the stars are taken 
from UV data. The open triangles are from PBSTG and are in common with our sample. 



Table 1. Observed Planetary Nebulae 
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a Planetary nebulae morphology: R = Round, E = Elliptical, and B = Bipolar. 
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Table 2. Observing Log 







Frequency 
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Av 


Telescope 


Transition 


(GHz) 


(arcsec) 


(km sec -1 ) 
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Table 3. Results 



Source 


Transition T peak 
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0.12 
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■*1 
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0.02 
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+1 


0.04 


0.01 
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^2 
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0.13 


i3CO J=3- 


^2 
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0.08 



12 CO J=2- 
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12 CO J=2- 
13 CO J=2- 



0.17 
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0.01 

0.05 
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12 CO J=2^1 0.12 0.01 
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-19.5 
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25.5 


8.8 ± 0.2 
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23. ± 0.9 


35.1 
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1.1 ± 0.1 
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10. l b 


18. ± 0.2 


31. 8 b 
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9.8 ± 0.1 


26.2 


49. ± 0.8 


31.1 


21. ± 1.0 


3.25 


0.92 ± 0.1 


10.7 b 


18. ± 0.4 


5.16 b 


2.6 ± 0.1 
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Table 3 — Continued 



Source 


Transition 


T 

- 1 peak 


RMS 


^centroid 




J 






(K) 


(K) 


( km sec -1 ) 


( km sec -1 ) 


(K km sec -1 ) 



NGC 6720 13 CO 


J= 


=2- 


■+1 


0.05 


0.01 


-3.5 


49.1 


1.2 


± 


0.2 


12 CO 


J= 


=3- 


^2 


1.4 


0.11 


-4.7 


16.3 


22. 


± 


0.9 


13 CO 


J= 


=3- 


^2 


<0.12 


0.04 













12 CO 


J=2- 


+1 


8.4 


0.003 


26.8 


22.3 


200. 


± 2.0 


NGC 7027 13 CO 


J=2- 


+1 


0.37 


0.004 


22.8 


27.4 


6.4 


± 0.9 


12 CO 


J=3- 


^2 


13. 


0.65 


26. 


22. 


300. 


± 2.0 


13 CO 


J=3- 


^2 


<1.0 


0.34 











12 CO 


J=2- 


+1 


2.5 


0.06 


18.1 


18.4 


50. ± 0.4 


IRAS 21282 13 CO 


J=2- 


+1 


0.06 


0.01 


15.6 


18.8 


0.66 ± 0.1 


12 CO 


J=3- 


^2 


5.1 


0.13 


18.9 


17.3 


100. ± 0.9 


13 CO 


J=3- 


^2 


0.15 


0.04 


21.8 


19.8 


1.5 ± 0.3 



12 CO J=2- 


+1 


3.3 


0.05 


-15.4 


3.5 


15. ± 0.2 


13 CO J=2- 


+1 


0.34 


0.02 


-15.3 


3.2 


1.6 ± 0.1 


i2 C J=3- 


^2 


0.97 


0.12 


-14.8 


2.6 


2.5 ± 0.2 


i3CO J=3- 


^2 


0.29 


0.08 


-15.7 


2.3 


0.52 ± 0.2 



a Based on a one-component Gaussian model. 

b Multiple components exists within the single Gaussian fit. 
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Tablc 4. 12 C/ 13 C Abundance Ratios 



Source 


PK 




12Q/13 


C Ratio 




NRAO 12 m 
(J=2^1) a 


HHT 10 m 

(J=3^2) a 


Adopted b 


PBSTG C 


AFGL 618 


(166-06) 


6.8 (6.9) 


5.2 (5.5) 


> 4.6 d 


4 e 


Ml-7 


(189-07) 


20. (18.) 


>6.7 


20. ± 1.8 




NGC 2346 


(215+03) 


22. (28.) 


>6.3 


22. ± 2.7 


23 


Ml-16 


(226+05) 


2.6 (2.7) 


2.4 (2.2) 


2.2 ± 0.03 d 


3 e 


VV 47 


(164+31) 


> 5.7 




> 5.7 




NGC 6072 


(342+10) 


6.9 (4.2) 




6.9 ± 0.31 




M4-9 


(024+05) 


> 4.0 




> 4.0 


18 


NGC 6720 


(063+13) 


9.5 (10.8) 


> 14. 


9.5 ± 1.6 


22 


NGC 7027 


(084-03) 


31.(22.) 


> 13. 


31. ± 0.62 


25 c 


IRAS 21282 




76. (42.) 


66. (33.) 


> 32. d 




NGC 7293 


(036-57) 


9.4 (9.7) 


4.8 (3.3) 


9.4 ± 0.60 


9 



a Based on CO integrated intensities and peak intensity (in parentheses). 
b Adopted 12 C/ 13 C ratio (see text). 
c From Palla et al. (2000). 

d Based on LVG model. 

e Lower limit based on optically thick CO. 



